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We discuss an astrometric timing effect on data analysis of continuous gravitational waves from 
rapidly rotating isolated neutron stars. Special attention is directed to the possibility of determin- 
ing their distances by measuring the curvature of the wave fronts. We predict that if continuous 
gravitational waves from an unknown neutron star with a stable rotation are detected around 1kHz 
within l/3yr by initial LIGO detectors and the ellipticity parameter e is smaller than 10 -6 , the 
distance r to the source can be estimated with relative error Ar/r of ~ 10% by using the broad 
band configuration of advanced LIGO detectors over 3 years. By combining the observed amplitude 
l^-^ ■ of the waves with the estimated distance, information on the parameter e can be obtained purely 

' through gravitational wave measurements. 
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I. INTRODUCTION 

Over the last 5 years, several large-scale interferometric gravitational wave detectors, US LIGO Q, Japanese TAMA 
, and British-German GEO , were constructed and have given us scientific data. Italian- French VIRGO £| will be 
soon in operation. In the next decade, second generation detectors (e.g. advanced LIGO @, LCGT Q) are expected 
to be available. More ambitious third generation detectors (e.g. EURO 0) will also be realized in the future. While 
iy-j ■ the current central objective of scientific runs is the direct detection of gravitational waves in 10Hz to a few kHz 
■^j- | band, gravitational wave astronomy will be rapidly established with the advent of new and more sensitive detectors 
i 0,0- Once gravitational waves from a population of sources are detected with relatively small signal-to- noise ratios 
(SNRs), follow-on detectors would reveal new astronomical information at high SNRs. This kind of consideration is 
important for the design and operational strategy in planing next generation detectors. 

In this paper, we study an astrometric effect relevant to the data analysis of gravitational waves from rapidly 
rotating isolated neutron stars and discuss the possibility of estimating their distances. If gravitational waves from 
■ a rotating neutron star are detected, the amplitude contains information about the distance r and the ellipticity 
Oh | parameter e through the combination e/r. The parameter e characterizes the non-axisymmetry of the star. However, 
it is not easy to separate these two quantities individually from gravitational waves. Here, we show that by a long 
term observation of gravitational waves from a source, we can measure a small phase shift induced by an astrometric 
effect. This means that a rapidly rotating neutron stars are the ideal targets to apply the long-term astrometric 
studies with planned ground based detectors. Gravitational waves from merging neutron stars or supernovae would 
\ last at most order of seconds 0,0 and are not suitable for detecting the shift. It is also interesting to understand how 
the quality of the astronomical information depends on the observational duration. For a very short phenomenon, 
the only information improved by a long data span is the event rate. 

To confirm gravitational waves from a rotating neutron star, we have to integrate the signal for a long time to 
get a sufficient signal-to- noise ratio [ToL [TT| (see 01 f° r a recent observational study with GEO and LIGO). The 
required computational cost for the data analysis depends strongly on whether the neutron star is known or unknown. 
If its location is unknown, we need a large number of templates, and the data analysis is a very challenging task 
[l3L Il4|. The astrometric effect is larger for a closer object, and we implicitly assume a nearby unknown neutron 
star as our main target. Once gravitational waves from an unknown neutron star are detected, we might identify an 
electromagnetic- wave counterpart using information obtained from the gravitational wave measurement, and estimate 
its distance with various traditional methods (see e.g. |l5fV But a search for an electromagnetic- wave counterpart 
would be very difficult for a cold neutron star with a large peculiar velocity and a very weak magnetic field (e.g. less 
than ~ 10 7 gauss critical for millisecond pulsars [T^h 

This paper is organized as follows: In the next section we formulate the phase shift induced by the astrometric 
effect for a rapidly rotating neutron star, and explain a way to estimate its distance though the phase shift. Then, 
using the Fisher matrix approach, we obtain an expression for the magnitude of error associated with this estimation. 
This expression is our central result. In section III we apply it to observational situations with specifications for future 
projects, such as, advanced LIGO. We examine required conditions for the ellipticity e and the source distance r in 
order to detect gravitational waves from rotating neutron stars and to determine their distances with the proposed 
method. Then we statistically study whether there are nearby neutron stars to apply our method. Section IV is 
devoted to a discussion including a brief summary. 



2 



II. FORMULATION 



We first study how well we can extract information on the distance from the phase modulation of gravitational waves. 
Similar arguments have been made for radio pulsar timing analyses [l7l Il9j , but this is the first attempt at applying it 
to gravitational wave data. Our basic aim is to develop a model relevant for evaluating parameter estimation errors. 
We do not try to make a detailed model required for actual signal analysis e.g. including general relativistic time 
delays in the solar system that are small and introduce no additional source parameters (Tj. In addition, the effects 
of diurnal rotation are not included here, because we are dealing with the case in which the signal integration is of 
order a few years. Following Ref.|ll|], we use the solar system barycenter (SSB) rest frame and investigate the effects 
induced by the motions of the detector as well as the source. We denote the time dependent position of the detector 
I'd {rd = \rd\ ~ 1AU), and assume that the source is moving with a constant velocity v ns as rn,Q + v ns (t + r/c) 
in the SSB frame. Here too is the position of the source (no; source direction) from which the emitted wave front 
reaches the SSB at t — 0. We define no = (sin 9 cos <j>, sin 9 sin cf>, cos 9) in the spherical SSB coordinate with 9 being 
the angle between the source direction no and the normal vector of the orbital plane (ecliptic) of the Earth. We write 
the gravitational wave phase ^ssb observed at the SSB in the form ^sSB(t) = $o + 2-7T J2k=o fkt k+1 /{k + 1)! and 
perturbatively include the phase modulation due to the detector's and source's motions with two small expansion 
parameters; a — \rd/r\ and (3 = \v ns t\/r. After some algebra, the phase ^d observed by the detector is given to 
0(a m (3 n ) {n + m < 2) as 
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k=0 y ' v 7 fe=0 

+higher order terms, (2-1) 

where v ns ± is the velocity component of v ns perpendicular to the direction no |l7j . There are no terms of order 
0((3) and 0{(i 2 ) EJ. In the large parenthesis in the above equation, the O(a) term is nothing but the plane wave 
effect that has been analyzed in the literature and can be used to determine the direction of the source. This term is 
periodic with a lyr -1 frequency. The 0(a(3) term was discussed to some extent for gravitational wave data analysis 
in Ref. [TTL |20| and provides information on the proper motion of the source on the sky. The distance r to the source is 
estimated through the 0(a 2 ) term. This term probes the curvature of the wave fronts as they deviate from plane wave 
propagation. In the 0(a 2 ) term, the constant contribution — l/(2r 2 ) is effectively absorbed in the phase constant $0, 
while the time dependent oscillating term (with a 2yr _1 frequency) is proportional to sin 2 9. The order of magnitude 
of the time shift induced by the plane wave 0(a) term in one year is written as ~ r d cos0/c ~ 500sec with G being 
the angle between the two vectors no and rj,. This angle changes by AO ~ rd/r in one year due to the parallax. 
Therefore the 0(a 2 ) term ~ r 2 d /cr ~ 10 -4 (r/30pc) -1 sec can be regarded as a time shift induced by the parallax. Even 
at / ~ 1kHz, this time shift produces a phase shift much smaller than the 1/4 cycle that is critical for detection of 
gravitational waves, unless the source is very close to the solar system. Therefore the distance r cannot be estimated 
without a high SNR observation well beyond the required level for gravitational wave detection. 

Using a model h — h c sm^d] (h c - constant) for the gravitational wave signal, we calculated the Fisher matrix for 
the following 10 fitting parameters {fo, /ij jii /3: ^o> Vns,±, r i $7 and examined the relative estimation error Ar/r 
with various sets of input parameters as well as the total observational time T Q b s > lyr. We found that the result Ar/r 
is independent of the parameters 4>, $0 and inputs {/1, /2, /3, /4, v ns ±}, as long as these inputs are close to reasonable 
values for real neutron stars. For a integration time T a b s longer than ~ 2yr, the resolution Ar/r is approximated well 
by the following expression 

Ar . / r \ r SNR y 1 ( sin6> \~ 2 / / 



' U VIOOpcJ \ 500 ) Vsin7r/3j ^lkHzJ ' ^'^ 

We take 9 — ir/3 for a reference value for the angle 9, as it bisects the area of a half sphere < 9 < tt/2. With the 
specific choice of parameters 9 — tt/3, f = 1kHz, r = lOOpc and SNR — 500, the prefactor 0.11 in eq.(2) becomes 
9.14, 0.13, 0.109, and 0.107 for T f, s = 1,2,3 and 4 yr respectively. This indicates that for T b s > 2yr the correlation 
between r and other parameters decreases significantly and becomes nearly stationary. The coefficient becomes 0.107, 
when we remove f% or v ns ± from the fitting parameters with a fixed observation period T f, s = 3yr. In our model 
we have not directly dealt with the effects of the acceleration of the source. However, these numerical experiments 
imply that, even if we include acceleration in a relatively simple manner, it would not change the resolution Ar/r 
significantly. When the source is a binary (e.g. with its orbital period ~ lyr), the situation could change considerably. 
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III. PROSPECTS WITH PROPOSED DETECTORS 

Next, we discuss how well we can estimate the distance to a rotating neutron star with the planned gravitational 
wave detectors, such as, advanced LIGO. We first describe the standard detection criteria of its gravitational wave 
signal based on SNR for a matched filtering method 0, ^| . Then we use a expression for SNR to evaluate the 
distance estimation error given by ea. l|2.2(l . For these studies, we pay a close attention to dependence of the ellipticity 
parameter e and the distance r to the source. 

If the rotation axis of the neutron star is identical to one of the principal axes of the inertial moment tensor , 
the orientation averaged amplitude of gravitational waves becomes h c — 8ir 2 Gy / 2/15If 2 e/c i r, where gravitational 
ellipticity is defined as e = {I22 — hi)/! <C 1, and / = ^33 (2:3: the rotation axis) |g. With typical numerical values, 
we obtain 




Hereafter we fix the inertial moment at / = 10 45 g cm 2 and the gravitational wave frequency at 1kHz. When the 
rotation axis is different from a principal axis, the star precesses and emits gravitational waves at multiple frequencies 
0, . In this paper we use the above equation (|3.1|1 as a reference to relate the characteristic amplitude h c with 
two parameters e and r. For the sensitivity of initial and advanced LIGO at 1kHz, we adopt the numerical values 
given in Figure 1 in Ref.@ (see also [3). We take h d = 1.5 x lO^H?" 1 / 2 for a single 4km detector of initial LIGO 
and hd = 4.2 x 10 _24 Hz -1 / 2 for advanced LIGO with a broad band configuration. After averaging with respect to 

the direction of the source, the effective sensitivity for a monochromatic signal becomes h e ff — h^yj 2~ 1 5T~ b ] for an 

observational period T b s . The factor 2 -1 / 2 comes from the number of 4km detectors and another factor 5 1 / 2 is due 
to the angular average with respect to the source direction. We also study the case for the proposed xylophone-type 
EURO detector Q with hd = 3.6 x 10 _26 Hz -1 ^ 2 at 1kHz. With these concrete specifications, the signal to noise ratio 
is calculated as 

SNR=-^, (3.2) 

h eff 

and the resolution Ar/r is now given as a function of r and e from eas. <|2.2[)i|3.1|) and l|3.2|) . 

The characteristic SNR threshold required for detecting an unknown rotating neutron star is discussed in |10| and we 
take SNR = 20 as a reference value for a l/3yr integration. In figure 1 we plot the r — e relation for four observational 
conditions; SNR = 20 with (i) initial LIGO and (ii) advanced LIGO for l/3yr integrations, and Ar/r = 0.1 with (iii) 
advanced LIGO and (iv) xylophone-type EURO for 3yr integrations. We have r oc e for a given SNR, but r oc e 1 / 2 
for a given resolution Ar/r. If the parameter e is fixed, we have SNR oc r _1 and Ar/r oc r 2 . 

To detect gravitational waves from neutron stars, a longer integration time is statistically advantageous. This 
is because of the fact that the amplitude of the detectable signal decreases and the effective survey volume (or 
equivalently the event rate) increases. However, the number of templates required for the matched filtering method 
depends strongly on the integration time, and a manageable data length is limited by computational resources. In 
reality, it might be difficult to deal with data of 1/3 years, especially for a young neutron star with large coefficients 
fk (k > 1) (see ea. ^.ll) 1 ). In this case, we can detect neutron stars that are closer than the distances shown by the 
dotted lines in figure 1. However, with a total ~ 10 Teraflops computational power and efficient searching algorithm, 
the differences between the dotted lines and the maximum detectable distances for initial and advance LIGO would 
be less than a factor of 2 for most neutron stars (see figure 5 in ^3)- Once gravitational wave from a neutron 
star is confirmed, it would be a relatively easy task to analyze longer data streams with a finer spacing of template 
parameters. But the available data span would not be much larger than 3 years, considering the time scale of major 
upgrades to next generation detectors. This is the reason why we studied the case with a 3yr integration for a follow-up 
analysis of the distance measurement. 

In figure 1, the solid line for advanced LIGO and the dotted line for initial LIGO intersect around the interesting 
value e = 10 -6 , close to the upper limit from theoretical predictions |l0j. Almost all the triangle region in (e, r)-plane 
bounded by e < 10~ 6 (dashed line) and the dotted line for initial LIGO is below the sold line for advanced LIGO. 
From viewpoints of future prospects of the gravitational wave astronomy, this fact can be rephrased as follows: if an 
unknown neutron star is detected around 1kHz with initial LIGO by ~ 3 _1 yr signal integration and the ellipticity has 
a reasonable magnitude e < 10 -6 , advanced LIGO can determine its distance with error of ~ 10% level within 3 years. 
This is quite encouraging because blind searches for neutron stars are very hard tasks requiring a huge computational 
cost, as mentioned earlier Once it is detected with initial LIGO, advanced LIGO would provide information 

about its distance, one of the fundamental astronomical parameters for a source. The ellipticity parameter e may 
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FIG. 1: Relation between the distance r and the gravitational ellipticity e for given observational conditions of an unknown 
neutron star with gravitational wave frequency / = 1kHz. The thin dashed lines represent the maximum distance for detection 
of gravitational waves with SNR=20 using initial LIGO and advanced LIGO for l/3yr integrations. The thick solid lines 
are the maximum distance for resolution Ar/r = 0.1 through the timing parallax method, shown for advanced LIGO and 
xylophone-type EURO for 3yr integrations. The dashed line e = 10 -6 is an upper limit of theoretical predictions for the 
ellipticity parameter. 



then also be found through knowledge of the amplitude h c . Both parameters are thus determined purely through 
gravitational wave measurements. Note that advanced LIGO has the potential to measure the distance to an unknown 
neutron star that cannot be detected by initial LIGO. The effective sensitivity of the single EURO xylophone- type 
detector is 82 times better than advanced LIGO detector at 1kHz, and the distance r for a given resolution Ar/r is 
9.1 times larger. 

We have discussed the case in which the gravitational wave data would be integrated without a break (100% duty 
cycle) for 3 years or so. But our result for Ar/r would be approximately valid for a coherent analysis of intermittent 
data streams that are sampled with appropriate intervals (<C lyr) and in a sufficient span (e.g. 3 years). If the narrow 
band configuration is available at the target frequency of a source, we can reach the same resolution Ar/r with a 
shorter total integration time. For example, the sensitivity hd could be improved by a factor ~ 3 at 1kHz with the 
narrow band configuration of advanced LIGO, and the total integration time could be a factor ~ 10 times smaller 
than for its broad band configuration. 

From Figure 1, we can estimate the distance r for the required resolution Ar/r ~ 0.1 as a function of the parameter 
e. Here let us statistically study whether there is at least one rapidly rotating pulsar within this distance. For 
this purpose we use the argument originally made by Blandford assuming a population of neutron stars whose spin 
evolutions are primarily determined by energy loss through gravitational radiation [g, llOj - We define t^ 1 as the 
formation rate of the population in our Galaxy with an initial gravitational wave frequency higher than 1kHz. The 
magnitude of the time tb is highly uncertain but should be much larger than tb ~ lOOyr corresponding to the Galactic 
supernova rate. We first evaluate the minimum distance r m i n to a pulsar of this population, and then compare r m i n 
with the distance allowed for the resolution Ar/r = 0.1. If these two distances coincide, we can expect to have at 
least one pulsar with resolution Ar/r better than 0.1. Our goal here is to calculate the required formation rate tb 
for this condition as a function of the ellipticity e. 

In Blandford's model, the amplitude h c from a neutron star with a fixed distance becomes larger for a larger e. On 
the other hand, the evolution time scale tgw becomes smaller for a larger e, and the smallest distance r m i n would be 
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larger. Here the evolution time scale tqw is defined as 

r ™ - r = L8 x 108 (up*) " 2 ( i4) 4 yr - (3 - 3) 

The gravitational frequency becomes smaller by ~ 30% within a single tgw- The time scale tqw becomes the age 
of the universe ~ 10 10 yr at e ~ 10~ 9 (//lkHz)~ 2 . For the spatial distribution of the population in our Galaxy, we 
assume a uniform density in a cylinder with radius R ns — lOkpc and height H ns = lkpc. If the parameter e is in the 
range 10" 9 (//lkHz) 2 < e < 10^ 6 (r B /200yr)- 1 / 2 (//lkHz)^ 2 , we have 

" <**»*W " S " 44 fe) US (kr) 2 ' 3 (use) » <"> 

for the smallest distance. The distance r m i n becomes larger than the height H ns for e > 
10 _6 (r B /200yr)^ 1 / 2 (//lkHz)- 2 . In this case we have r m i n ~ -RnsfTBTcw) ^ £l ^ 2 - With these expressions and 
results presented in figure 1, we can now estimate the required formation rate tb for observing at least one pulsar 
with resolution Ar/r = 0.1. For 3yr advanced LIGO observation at 1kHz, we have t# < 84(e/10 _8 ) _1 / 2 yr in the 
range 10~ 9 < e < 10~ 5 . Therefore, a very small tb (or a large formation rate t^ 1 ) is required from this rough 
statistical argument. In the case of 3yr EURO-xylophone-type data, we obtain t b < 6.3 x lO^e/lO" 8 )" 1 ^ for 
10~ 9 < e < 10- 7 . 

So far, we have studied the case with Ar/r <C 1. Here, we discuss the opposite case in which the distance to a 
source is large and we cannot detect the signature of the parallax timing shift. In the case of usual angular parallax 
measurement, the amplitude of the apparent annual motion of a source on the sky is inversely proportional to its 
distance r. Therefore, it is straightforward to fit distance with a form r _1 in the parameter estimation. We can 
understand this from the original meaning of the distance unit "parsec" that represents 1/arcsec. The situation is 
similar to the fitting for the parallax timing. Note also that the error A(r _1 ) ~ (Ar)/r 2 for r _1 does not depend on 
the distance r itself, as shown in eq.(2). In the low resolution limit with a fitting result |(?" -1 )/it| ~ A(r _1 ), we can 
set a constraint for the value r -1 in the form r" 1 < max{(r _1 )/i t , 0} + A r A(r~ 1 ) with N — a statistical significance. 
At the end we get a lower limit of the distance r as 

for T b a > 2yr with x = max{(r^ 1 )/i t /A(r^" 1 ), 0}. We can also set a lower limit to the parameter e by using the 
gravitational wave amplitude h c . 



IV. DISCUSSIONS 



In our analysis, the phase $sss(i) at the SSB is assumed to be simple and be described by a low-order Taylor 
expansion. Even though many pulsars are known to be extremely accurate clocks, the phase $sss could have a 
deviation from this simple expression. The origin of these timing deviations can be due to the intrinsic evolution of 
the pulsars or other effects {e.g. the low frequency gravitational wave background). If the deviation is large, especially 
around the frequency ~ 2yr _1 , it would hamper the estimation of the distance through parallax timing even with a 
high SNR measurement. Here, we briefly describe the actual radio observation of certain pulsars following [2^]. PSR 
B1855+09 is a binary millisecond pulsar with spin period 5.36 msec and orbital period 12.3 day, and PSR B1937+21 
is a single millisecond pulsar with spin period 1.56 msec. They are considered to be recycled pulsars. Timing analyses 
for both pulsars have been performed with data taken at Arecibo Observatory in a ~ 8 year span. A systematic timing 
deviation of order a few /xsec is indicated for PSR B1937+21, but not for PSR B1855+09. The timing residual of the 
former is expected to be dominated by the intrinsic evolution of the pulsar and not by other effects, such as, the low 
frequency gravitational wave background. As for the estimated distances through the timing parallax, PSR B1855+09 
has a fitted value tt = (lkpc/r)mas = 1.1 ± 0.3mas corresponding to r = O^Q^kpc (lcr), and PSR B1937+21 has a 
bound 7r < 0.28mas (2a) corresponding to r > 3.6kpc. Their distances are also estimated from dispersion measures 
together with a model of the Galactic distribution of interstellar free electrons. They are 0.7kpc for PSR B1855+09 
and 3.6kpc for PSR B1937+21 with uncertainties of ~ 25%, and are consistent with the estimation through the timing 
parallax. Therefore, the effects of the timing deviation would be moderate in these cases. 

Now we summarize this paper. We have discussed an astrometric effect for the analysis of gravitational waves that 
would allow us to determine the distance to a rapidly rotating isolated neutron star around / ~ 1kHz. The resolution 
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Ar/r improves significantly, if the integration period of the signal is longer than ~ 2 years. By using advanced LIGO 
detectors for 3 years at its designed sensitivity in the broad band configuration, we can estimate its distance with an 
error of Ar/r < 0.1 up to r ~ 150(e/10" 7 ) 1 / 2 pc, if its spin evolution is relatively simple. 

Acknowledgments: The author is grateful to Teviet Creighton for valuable discussions. He also thanks Sterl Phinney, 
Kip Thorne and Yi Pan for their help and useful conversations, and Asantha Cooray and Steve Furlanetto for carefully 
reading the manuscript. This work was supported by NASA grant NNG04GK98G and the Japan Society for the 
Promotion of Science. 



A. Abramovici et al., Science 256, 325 (1992). 

M. Ando et al. [TAMA Collaboration], Phys. Rev. Lett. 86, 3950 (2001) |arXiv:astro-ph/0105473l . 

B. Willke et al, Class. Quant. Grav. 19, 1377 (2002). 

A. Giazotto, Nucl. Instrum. Meth. A 289, 518 (1990). 

E. Gustafson, D. Shoemaker, K. Strain and R. Weiss, 1999, LSC White Paper on Detector Research and Development 

(LIGO project document T990080-00-D) 

K. Kuroda et al, Class. Quant. Grav. 20, S871 (2003). 

EURO homepage: http:/www. astro. cf.ac.uk/geo/euro 

K. S. Thorne, in Three hundred years of gravitation, edited by S. W. Hawking and W. Israel (Cambridge University Press, 
Cambridge, 1987), pp. 330-458. 

B. F. Schutz, Class. Quant. Grav. 16, A131 (1999) arXiv:gr-qc/9911034 . 

P. R. Brady, T. Creighton, C. Cutler and B. F. Schutz, Phys. Rev. D 57, 21 01 (1998) |arXiv:gr-qc/9702050| . 
P. Jaranowski, A. Krolak and B. F. Schutz, Phys. Rev. D 58, 063001 (1998) arXiv:gr-qc/ 98040 14 . 
B. Abbott et al. [LIGO Scientific Collaboration], Phys. Rev. D 69, 082004 (2004) arXiv:gr-qc/0308050 . 
P. R. Brady and T. Creighton, Phys. Rev. D 61, 082001 (2000 ) |arXiv:gr-qc/9812014| . 

P. J aranowski and A. Krol ak, Phys. Rev. D 61, 062001 (2000) arXiv:gr-qc/9901013 ; C. Cutler, I. Gholami and B. Krish- 
nan, arXiv:gr-qc/0505082 

D. A. Frail and J. M. Weisberg, Astron. J. 100, 743 (1990); G. C. Gomez, R. A. Benjamin, D. P. Cox and P. Donald, 
Astron. J. 122, 908 (2001). 

E. S. Phinney and S. R. Kulkarni, Annu. Rev. Astron. Astrophys. 32, 591 (1994). 

R. Blandford, R. W. Romani and R. Narayan, Journal of Astrophysics and Astronomy, 5, 369 (1984); T. Damour and N. 
Deruelle, Ann. Inst. Henri Poincare 44, 263 (1986); D. C. Becker and R. W. Hellings, Annu. Rev. Astron. Astrophys. 24, 

537 (1986). 

A. Buonanno and Y. b. Chen, Phys. Rev. D 64, 042006 (2001) |arXiv:gr-qc/0102012| . 

L. A. Rawley, J. H. Taylor and M. M. Davis, Astrophys. J. 326, 947 (1988); M. F. Ryba and J. H. Taylor, Astrophys. 
J. 371, 739 (1991). 

P. Jaranowski and A. Krolak, Phys. Rev. D 59, 063003 (1999) arXiv:gr-qc/9809046 . 

M. Zimmermann and E. Szedenits, Phys. Rev. D 20, 351 (1979) - 

V. M. Kaspi, J. H. Taylor and M. F. Ryba, Astrophys. J. 428, 713 (1994). 



